An Effective Model of Magnetoelectricity in Multiferroics RMn^O^ 
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An effective model is developed to explain the phase diagram and the mechanism of magnetoelec- 
tric coupling in multiferroics, RMnzOs. We show that the nature of magnetoelectric coupling in 
RMniOh is a coupling between two Ising-type orders, namely, the ferroelectric order in the b axis, 
and the coupled magnetic order between two frustrated antiferromagnetic chains. The frustrated 
magnetic structure drives the system to a commensurate-incommensurate phase transition, which 
can be understood as a competition between a collinear or col-plane order stemming from the 'or- 
der by disorder' mechanism and a chiral symmetry order. The low energy excitation is calculated 
and the effect of the external magnetic field is analyzed. Distinct features in the electromagnon 
spectrums in the incommensurate phase are predicted. 

PACS numbers: 77.80.+q, 75.47.Lx, 77.80.-o 
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Recently, the search for new spin-electronics materials 
has led to a discovery of novel gigantic magnetoelectric 
and magnetocapacitive effects in rare-earth manganites, 
magnetoelectric multiferroics [3, Hj]. Unlike the mag- 
netic ferroelectroics studied in 1960s and 1970s where 
magnetism and ferroelectricity couple weakly, the mag- 
netism and ferroelectricity in the new materials couple 
so strongly that the ferroelectricity can be easily ma- 
nipulated by applying a magnetic field and the magnetic 
phase can be controlled by applying an electric field [H, 0|. 
This ease of manipulation promises great potential for 
important technological applications in novel spintronics 
devices. 

The physics of the multiferroics involves the inter- 
play between many degrees of freedom, such as charge, 
spin, orbital and lattice. Tremendous effort has been 
devoted to decode the fundamental mechanism of the 
strong coupling between the magnetism and ferroelec- 
tricity. Experimentally, two major classes of magnesium 
oxide multiferroics, have been discovered. The first class 
is the orthorhombic rare-earth manganites RMnO^iR — 
Gd,Tb,Dy, ...)@, H[, characterized by spiral magnetism 
strongly coupled with the ferroelectricity. An effective 
Ginzburg-Landau theory incorporating the space group 
symmetry and time reversal symmetry has been con- 
structed to explain the fundamental physics [7J. Micro- 
scopically, Dzyaloshinskii-Moriya spin-orbit interaction is 
the underlying mechanism of the ferroelectricity [^.lilficij]. 
The second class of materials are the manganese ox- 
ides with general formula RMn20^{R — Y,Tb, Dy, ...) 



Ill |12| . Il3l . Il4| . These insulating materials consist of 
linked Mn i+ C>Q octahedra and Mn 3+ Os pyramids with 
a Pbam space group symmetry. Unlike that in RMnO^, 
the ferroelectricity in RMriiOe, exists in a collinear or col- 
plane magnetic phase, suggesting that a different mecha- 
nism is involved in the interaction between the ferroelec- 
tricity and magnetism. 

In this Letter, we develop an effective model to explain 
the phase diagram and the mechanism of magnetoelec- 



tric coupling in RMn-iO^. Building upon experimental 
facts and the space group symmetry [J, [ll|, [l2], |l3j, llJ, 
SHE EE Hi, we show that the magnetoelectric 
interaction is between two Ising type orders, the ferro- 
electric order in the b axis and the coupled magnetic 
order between two frustrated antiferromagnetic chains. 
The effective model of the magnetism can be derived 
from a microscopic model with nearest-neighbor mag- 
netic exchange. We show that the effective model nicely 
captures the phase diagrams of RM^O^. At high tem- 
perature, the commensurate (CM) collinear or col-plane 
order is stable due to the 'order by disorder' mechanism 
(2fl 21, 22 1 and the existence of an easy axis. As the 



temperature decreases, a chiral symmetry order replaces 
the collinear or col-plane order, and the magnetic struc- 
ture becomes incommensurate (ICM). This model pre- 
dicts that an external magnetic field along the b-axis can 
drive the system from the ICM phase to the CM phase 
while that along the a-axis can drive the system from 
the CM phase to the ICM phase. The model predicts 
the following distinguished properties of low energy exci- 
tations: (1) the emergence of electromagnons in the ICM 
phase when there is no electromagnons at the lowest en- 
ergy in the CM phase; (2) the presence of distinguished 
kinks in the energy dispersions of the electromagnons, 
unlike the dispersion of normal phason in conventional 
ICM phase which has a finite energy jump at the half 
of ICM wavevector [23[; (3) a double peak structure at 
low energy in optical conductivity due to the absorption 
of the electromagnons with a selection rule of the elec- 
tric field along the b axis; (4) no new peak splitting in 
the electromagnon spectrums in the presence of external 
magnetic field along the a or b axis, which differs from 
the conventional picture of the Zeeman energy splitting 
of magnons; (5) the increase (decrease) of the energy gaps 
of the electromagnons as the field increases along the b(a) 
axis. 

Magnetoelectric coupling: The ferroelectricity in 
RMn^O^ is substantially different from that in RMnOj,. 
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Experiments have shown that the ferroelectricity only ex- 
ists along the b axis in RMriiO^, but can be observed 
in both the a and c directions in RMnO^. Most im- 
portantly, recent measurements of optical conductivity 
have revealed opposite selection rules for these two ma- 
terials. In RMnOs, low energy absorption is observed 
when the electric field is perpendicular to the static fer- 
roelectric polarization direction ljj 24 1 . The opposite is 
true in RMn^O^^ namely, low energy absorption is only 
observed when the electric field is polarized along the b 
axis This suggests that the electric degree of free- 
doms along the a and c axis simply has no coupling to 
magnetic degrees of freedom at low energy. Only the 
ferroelectricity along the b axis Pb couples to the mag- 
netic degree of freedom. The order parameter Pb is an 
Ising-type order. Therefore, the nature of the magneto- 
electric coupling in RMTI2O5 is a coupling between two 
Ising-type orders. ( in the following paper, we refer the 
a,b,c axis to the x,y,z axis respectively for conveniences 



i.e. P v = Pb). 



What is the Ising order in the magnetic side? The an- 
swer to this question can be obtained by in-depth anal- 
ysis of the magnetic structure and the space group. As 
shown in 16|, 13], the main magnetic structure along the 
a axis is two antiferromagnetic chains joined by Mn 3+ 
and Mn 4+ atoms (see FigfTJ) . The antiferromagnetic cou- 
pling between the two chains indicated by the red lines in 
FigHJare completely frustrated. Therefore, in an effective 
model, we at least need two antiferromagnetic orders n\ 
and n,2 to describe the magnetic physics. A possible Ising 
order from these two vector magnetic orders is n\ ■ n^. 
Due to the experiment fact that no magnetic moment in 
the c-axis is observed, it naturally leads to the construc- 
tion of the possible lowest order magnetoelectric coupling 
1,2 as 



bewteen P y and n, , i 



V V 



(1) 



The possible differences between the coupling parameters 



A , 



and A y reflects real lattice structure. 



Now we show that EqQ] is consistent with the space 
group analysis. The space group of RMn 2 0^ has been 
analyzed [25[. The lattice of RM^O^ belongs to Pbam 
structure. With the modulation vector q = (1/2,0, k c ), 
the space group has a single two dimensional irre- 
ducible representation in which the four symmetry lat- 
tice transforms can be represented by I,m x — a x ,m y — 
a y ,m x m y — ia z where <Ji are Pauli Matrix. Symmetry 
adapted variables can be constructed as linear combi- 
nations of spin operators that transform in accordance 



with these matrices. The ion spins in one unit cell are 
numbered one to eight as shown in Fig[T] The space in- 
version symmetry, together with the experimental facts 
that Si — S3, S2 = Sa,S§ = Sj, Sq = —Ss and that 
the spin moments are only in a-b plane, suggests that 
the possibil ities of magnetoelectric coupling term can be 




FIG. 1: A sketch of spin structures from a top view along the 
c axis in one unit cell of RMniO^. The red lines reflect the 
frustrated magnetic coupling between two chains. 



narrowed down to 

H em = i\ x P y (-S%(q)SI*(q) + S?(q)S%*(q) ~ cc) 
+ iX v P y (-Sl( q )Sr(q)+Sl( q )Sr( q )-c.c). (2) 

Converting these spin operators to the two antiferromag- 
netic orders, we can simplify Eq|2]to EqfT] 

Effective magnetic model: In condensed matter 
physics, an effective model at low energy is largely in- 
dependent of microscopic models if they share same es- 
sential physics. Therefore, one can derive the effective 
model on a much simplified lattice structure. In the case 
of i?Mn 2 05, the important magnetic physics along the 
a axis are two antiferromagnetic chains with frustrated 
coupling[l^, 17 1 . Experiments have shown a 1/4 com- 
mensurate magnetic wavevector along the c axis, which 
can also be viewed as an antiferromagnetic order if the 
unit cell is doubled along the c-axis. Therefore, we can 
derive the effective model from a microscopic model with 
two antiferromagnetic orders defined on the two interpen- 
etrating sublattices as illustrated in Figj2l where Ji (2) 
are the effective antiferromagnetic exchange couplings 
which establish two antiferromagnetic orders, J3 is the 
effective frustrated coupling between two chains in one 
unit cell along the a axis and J4 is the effective frus- 
trated coupling between two chains in two neighbor unit 
cells along the c axis. Using standard field theory [22I |26j|. 
we can show that the effective field theory described by 
the two antiferromagnetic orders n± and n2 is given by 
the following Hamiltonian 
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2 -(dn l /dzY}+a(n 1 — 



dni 
n 2^r~ ) 



dm dn 2 dn x dn 2 . , , 2 

* ( fl^ + ^^ )_fl(T)(ni,B3) 



dx 
*) 2 }dxdz. 



(3) 



f- Jl --* 



• J: \ • • 



J2 



• p— • • 

• * 



FIG. 2: A sketch of lattice structures of two frustrated coupled 
antiferromagnetic Heisenberg models in two dimension. 



where g{T) = go + g\T is a temperature dependent 
parameter induced by the quantum and thermal fluc- 
tuation, the parameter a reflects that the intra frus- 
trated coupling in one unit cell is larger than the in- 
ter frustrated coupling along the a axis between two 
neighbor unit cells along the c axis, and the param- 
eter D describes a possible magnetic easy axis along 
the a axis. From the microscopic coupling parame- 

2A-, 



-II J'd + Ji \2 ~ 
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ter, g x = 0.26(Ji + J 2 )Sa~ 2 {^^ )*, g 2 - «.-* (Jl+Ja)s , 

n — 7 C"2 _ T C2 n _ '(Jl+J2)S 2 (J3-J4) „ n j 

pis — Jib , p 2s — J 2 b , a — 1 (j 1+ j 2 ) a ana 

v = {.h+.h)S 2 Adding the lattice dynamics, we 

reach the total effective Hamiltonian as 



2M^ y) 



H„ 



(4) 



where ir y is the conjugate momentum of P y . 

The Hamiltonian in EqfJ] precisely captures the phase 
diagrams of RMn 2 0$^, [ll|, 0|. To study the mag- 
netic phase diagram, we can integrate out the lattice 
dynamics. After integrating out the lattice dynamics, 
the effective magnetic Hamiltonian is the Hamiltonian in 
Eq 11 01 with a replacement of g by g(T) = go + g\T where 



5o = 9o + D by D = D 



A + A_ 
4k 
A 2 



and an additional 



term — ^(n\n 2 — riffi^) 2 where 7 = and \± = X x ±X y . 
It is clear that the a term favors an ICM phase while 
g(T),j and D favor a CM phase. If D ^ 0, at relatively 
high transition temperature T c ± , the model exhibits a 
first phase transition to a collinear magnetic phase with 
order < nfn 2 >^ and then exhibits a second phase 
transition at T c2 < T c \ with order < n\n 2 >+ i 0. The 
collinear phase becomes a col-plane phase. The ferroelec- 

tricity is given by < P y >= — — - — 3 — 2— . With 

the proper values of a, at a low temperature Tic, the 
ICM phase can win over the col-plane magnetic phase. 



In the ICM phase, the global average, < P y >— 0. Fig[3] 
sketches the phase diagram. The phase diagram quali- 
tatively matches the current experimental results on the 
phase diagram of RMn 2 0^ [1, [111, [TBI . Results from mean 
field or large N limit calculation, using the real experi- 
mental data as input, will be reported elsewhere [2 7j. In 
this paper, we present a thorough study of the simplified 
version of the model in one dimension which captures the 
essential low energy physics. 

Coupled Sine-Gordon model: We focus on an one- 
dimensional effective model by ignoring the dynamics in 
z direction in Eqf4] The one dimensional model still cap- 
tures essential physics of the frustration. For the one 
dimensional model, we can go beyond mean field calcula- 
tion and evaluate the dynamics in a controllable pertur- 
bation manner. Using the angle to parameterize the an- 
tiferromagnetic vector in x-y plane as Hi = (cos9i, sinOi), 
we obtain the effective Hamiltonian in one dimension as 
a coupled Sine-Gordon model, 



H 



pde_, 

2^ dx ' 
2 [ dx ■ 



a sin 0_ — 

ox 

D cos Qjl. cos 0_ 



g cos 



7 cos 



(5) 



where 9± = 9\ ± 9 2 . Minimizing the potential part in 
Eqf5l the CM-ICM phase transition line is estimated to 
be 



p(g(T) +1 /2 + D) = a 2 /2 + 



D 2 



2{c?-pg{T)y 



(6) 



To obtain a more accurate result, we can apply a 
controllable numerical method [jjl]. In the CM phase, 
the ground state of the model is always described by 
9 + = 8- = or 7T. In the ICM phase, it is clear 
that the spins on the two antiferromagnetic chains are 
both rotating with the same wave vector q through- 
out the system. Therefore, we can take the follow- 
ing ansatz: 6-(x) = X)^Lo(°™ C0Sn Q x + °n svanqx) and 
Q+(x) — qx + X)^Lo( c « cos nax + d n smnqx). The coef- 
ficient in this ansatz decreases rapidly as n increases. 
As an example, for a typical set of parameters {p = 
2, a = 0.4, ff = 0.03,7 = 0.01, d = 0.005}, a nu- 
merical solution of the ICM ground state is given by 
9+ = qx + 0.0322 sin(2 9 a;) - 0.000255 sin(4gx) + ... and 
9- = -§ +0.0526 cosmic) +0.0000676 cos(3ga;)-h.., where 
q = 0.199. Therefore, the first two terms in the ansatz 
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FIG. 3: The phase diagram as go vs T. The vertical dashed 
line represents phase transitions in most of KMn^O^, mate- 
rials. The inset (a) is the ferroelectricity as a function of 
temperature. 



can give us a very good approximation of the ground 
state. 

Energy dispersion of magnons in the CM phase: In the 
CM phase, the ferroelectricity, < P y >= — A+ 2^ A ~ ■ The 
dynamics of magnons and the dynamics of electric de- 
gree of freedoms are decoupled. To show this, one can 
expand the free energy at the vicinity of the ground state, 
u\ = 59 '_(_, u 2 = 59-, [3 — 6P y and show that there is no 
second order coupling between Ui and (3 (the lowest order 
coupling between m and (3 is \ + [3u 2 +X-f3uf). Therefore, 
there is no electromagnon excitation. The energy disper- 
sions of the two magnons in the CM phase are given by 



El 



\J p{pk 2 + 7 + .g + D± \/{g-"i) 2 + a 2 k 2 ). (7) 



Electromagnons in the ICM phase: In the ICM phase, 
< P y (x) >= — ^cos(qx) where q is the ICM wavector. 
The magnon and phonon dynamics are coupled due to 
the ICM modulation. Therefore, the low energy excita- 
tions are electromagnons. To calculate the dispersion, we 
expand the Hamiltonian at the vicinity of the ICM phase, 
ui = 59 + , U2 = 59- and (3 = 5P y . The fluctuation up to 
the second order of these dynamic variables is given by 



5H^ = SH { > + 5H\ >. The first term is the free part 
given by 

2 

fc 8=1 

with E 1 (k) = pk 2 +°§±+ A - (ai 4 A K +A - } , E 2 (k) = E 1 (k) + 

A 2 

Ao where Ao = aq — 2g — and the phonon frequency 
— \ffi (for convenience, we have taken the mass M = 1 
in Eq0]) . The second term is the interaction part 

5E[[ 2S> — J dx[aai cos(qx)u\u 2 + D sm(qx)uiu 2 



+~(3{—\-U\sin{qx) + X + u 2 )]. 



(9) 



In EqUJl the terms in the first line couple iti(fc) with 
u 2 (k± q) and vice versa and modify the gap Aq between 



two magnon modes. In general, as shown in FigHl the 
coupling creates a distinguished kink in dispersion curve 
around k — q/2 rather than a finite energy jump at 
k = q/2 which is normally expected in the ICM phase 
described in a simple Sine-Gordon model[23]. The terms 
in the second line describe the coupling between the 
magnon and the phonon. The coupling results in two 
general effects. First, the coupling leads to a change of 
phonon frequency. Under condition Wg >> pEi(0), the 
shift frequency of the phonon is roughly given by 5ui 2 = 

4(^-pB 2 (0)) + S(^-pE 2 (q)) - _ 

us to measure the magnons in optical conductivity. The 
optical conductivity is given by o~(tu) = ujhnlGppiu), 0)], 
where Gpp(ui,k) is the full propagator of the ft. In our 
model, we expect double peaks in the optical conductiv- 
ity at the gap energy of the two magnons in the ICM 
phase. Up to the second order, we have, Im[Gpp(ui, 0)] = 

7r[p§z5(u 2 P E 1 (q)) + p^5{u 2 



Second, the coupling allows 



pE 2 (0)) + (1 - 
In Fig (4) by numeri- 



cally solving the dynamical equations of 5H^ 2 \ we plot 
the result of cr(w) and the dispersions of electromagnons 
with a typical parameter setting {p, a, A+, A- , k, D, g} = 
{2, 0.4, 0.06, 0.06, 2, 0.002, 0.005}. 

Effects of external magnetic field: In our model, since 
n is the staggered moment field, the effect of an external 
magnetic field H is equivalent to creating an easy plane. 
However, since the staggered moment is in the a-b plane, 
the effect of the field is only important when the field is 
along the a or b axis. In the presence of the easy axis 
parameter D along the a axis, the effects of the external 
magnetic fields along the a axis, if a , and b axis, Hi,, have 
exact opposite effects. They simply change D to 



D(H a , H b ) = D- H 2 a /2p + Hl/2p 



(10) 



in EqU From E x (k), E 2 {k), Eq(9]and EqEl EqllO] leads 
to a few important and immediate predictions. First, the 
effects of H a and Ht do not depend on their directions 
along their own axis. Second, Hb can drive the system 
from the ICM phase to the CM phase while H a can drive 
the system from the CM phase to the ICM phase. This 
result has been observed experimentally in 12]. The crit- 
ical fields to drive the transition can be estimated for the 
one dimensional system from Eq|SJ Third, in the ICM 
phase, the energy dispersions of the electromagnons as a 
function of H a and iff, can be predicted. From E\(k), 
E 2 (k) and EqEl the energy of the electromagnons is ex- 
pected to increase (decrease) as Hb (H a ) increases. Fi- 
nally, the external magnetic field does not add additional 
peaks, which contradicts the conventional picture of the 
Zeeman energy splitting of magnons. 

In conclusion, we develop an effective model that ex- 
plains the phase diagram and the mechanism of mag- 
netoelectric coupling in multiferroics RMn 2 0§. To our 
knowledge, this is the first theoretical effective model 
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FIG. 4: The numerical result of the optical con- 
ductivity versus frequency with the parameters set to: 
{p,a, A+, X-,K,D,g} = {2,0.4,0.06,0.06,2,0.002,0.005}. 
The inset is the dispersion of the electromagons. 

for these materials. A detailed study of low energy ex- 
citations in one dimension is performed to explain the 
selection rules of electromagnon in optical conductivity 
measurements Our prediction of the electromagnon 
dispersion and its dependence on the external magnetic 
field can be tested in future experiments. A quantita- 
tive study of our model incorporating the experimental 
data on different RM^Os materials will be reported 
elsewhere [13]. We expect that the model presented here 
can be applied to other multiferronics materials where 
ferroelectricity is correlated to a collinear or col-plane 
magnetic phase. 
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Recently, the search for new spin-electronics materials 
has led to a discovery of novel gigantic magnetoelectric 
and magnetocapacitive effects in rare-earth manganites, 
magnetoelectric multiferroics [? ? ]. Unlike the mag- 
netic ferroelectroics studied in 1960s and 1970s where 
magnetism and ferroelectricity couple weakly, the mag- 
netism and ferroelectricity in the new materials couple 
so strongly that the ferroelectricity can be easily manip- 
ulated by applying a magnetic field and the magnetic 
phase can be controlled by applying an electric field [? ? 
] . This ease of manipulation promises great potential for 
important technological applications in novel spintronics 
devices. 

The physics of the multiferroics involves the inter- 
play between many degrees of freedom, such as charge, 
spin, orbital and lattice. Tremendous effort has been 
devoted to decode the fundamental mechanism of the 
strong coupling between the magnetism and ferroelec- 
tricity. Experimentally, two major classes of magnesium 
oxide multiferroics, have been discovered. The first class 
is the orthorhombic rare-earth manganites RMnO^{R — 
Gd,Tb, Dy, ...)[? ? ], characterized by spiral magnetism 
strongly coupled with the ferroelectricity. An effective 
Ginzburg-Landau theory incorporating the space group 
symmetry and time reversal symmetry has been con- 
structed to explain the fundamental physics [? ]. Micro- 
scopically, Dzyaloshinskii-Moriya spin-orbit interaction is 
the underlying mechanism of the ferroelectricity [111 
]. The second class of materials are the manganese ox- 
ides with general formula RMn20^(R — Y,Tb, Dy, ...) 
[? ? ? ? ] . These insulating materials consist of linked 
Mn i+ C>Q octahedra and Mn 3+ 0$ pyramids with a Pbam 
space group symmetry. Unlike that in RMnO^, the fer- 
roelectricity in RM112O5 exists in a collinear or col-plane 
magnetic phase, suggesting that a different mechanism 
is involved in the interaction between the ferroelectricity 
and magnetism. 

In this Letter, we develop an effective model to explain 
the phase diagram and the mechanism of magnetoelec- 



tric coupling in RM.niO*>. Building upon experimental 
facts and the space group symmetry [? ???????? 
? ], we show that the magnetoelectric interaction is be- 
tween two Ising type orders, the ferroelectric order in the 
b axis and the coupled magnetic order between two frus- 
trated antiferromagnetic chains. The effective model of 
the magnetism can be derived from a microscopic model 
with nearest-neighbor magnetic exchange. We show that 
the effective model nicely captures the phase diagrams 
of RMn^Os. At high temperature, the commensurate 
(CM) collinear or col-plane order is stable due to the 'or- 
der by disorder' mechanism [? ? ? ] and the existence of 
an easy axis. As the temperature decreases, a chiral sym- 
metry order replaces the collinear or col-plane order, and 
the magnetic structure becomes incommensurate (ICM). 
This model predicts that an external magnetic field along 
the b-axis can drive the system from the ICM phase to 
the CM phase while that along the a-axis can drive the 
system from the CM phase to the ICM phase. The model 
predicts the following distinguished properties of low en- 
ergy excitations: (1) the emergence of electromagnons 
in the ICM phase when there is no electromagnons at 
the lowest energy in the CM phase; (2) the presence of 
distinguished kinks in the energy dispersions of the elec- 
tromagnons, unlike the dispersion of normal phason in 
conventional ICM phase which has a finite energy jump 
at the half of ICM wavevector [? ]; (3) a double peak 
structure at low energy in optical conductivity due to the 
absorption of the electromagnons with a selection rule of 
the electric field along the b axis; (4) no new peak split- 
ting in the electromagnon spectrums in the presence of 
external magnetic field along the a or b axis, which dif- 
fers from the conventional picture of the Zeeman energy 
splitting of magnons; (5) the increase (decrease) of the 
energy gaps of the electromagnons as the field increases 
along the b(a) axis. 

Magnetoelectric coupling: The ferroelectricity in 
RMnyQ^ is substantially different from that in RMnO^. 
Experiments have shown that the ferroelectricity only ex- 
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ists along the b axis in RMn 2 0^, but can be observed 
in both the a and c directions in RAlnO^. Most im- 
portantly, recent measurements of optical conductivity 
have revealed opposite selection rules for these two ma- 
terials. In i?Mn0 3 , low energy absorption is observed 
when the electric field is perpendicular to the static fer- 
roelectric polarization direction [? ? ]. The opposite is 
true in RMnzOs, namely, low energy absorption is only 
observed when the electric field is polarized along the b 
axis[? ]. This suggests that the electric degree of free- 
doms along the a and c axis simply has no coupling to 
magnetic degrees of freedom at low energy. Only the 
ferroelectricity along the b axis Pb couples to the mag- 
netic degree of freedom. The order parameter Pf, is an 
Ising-type order. Therefore, the nature of the magneto- 
electric coupling in RM^O^ is a coupling between two 
Ising-type orders. ( in the following paper, we refer the 
a,b,c axis to the x,y,z axis respectively for conveniences 
i.e. P v = P b ). 

What is the Ising order in the magnetic side? The an- 
swer to this question can be obtained by in-depth anal- 
ysis of the magnetic structure and the space group. As 
shown in [? ? ],the main magnetic structure along the 
a axis is two antiferromagnetic chains joined by Mn 3+ 
and Mn 4+ atoms (see FigflJ . The antiferromagnetic cou- 
pling between the two chains indicated by the red lines in 
Figfflarc completely frustrated. Therefore, in an effective 
model, we at least need two antiferromagnetic orders n\ 
and n2 to describe the magnetic physics. A possible Ising 
order from these two vector magnetic orders is ni • n 2 . 
Due to the experiment fact that no magnetic moment in 
the c-axis is observed, it naturally leads to the construc- 
tion of the possible lowest order magnetoelectric coupling 
bewteen P y and fti , i = 1, 2 as 
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The possible differences between the coupling parameters 
X x and X y reflects real lattice structure. 

Now we show that EqQ] is consistent with the space 
group analysis. The space group of RMn 2 0$ has been 
analyzed [? ]. The lattice of RMn^O^, belongs to Pbam 
structure. With the modulation vector q — (1/2,0, k c ), 
the space group has a single two dimensional irre- 
ducible representation in which the four symmetry lat- 
tice transforms can be represented by I,m x — a x ,m y — 
<7 y ,m x m y — ia z where <Ji are Pauli Matrix. Symmetry 
adapted variables can be constructed as linear combi- 
nations of spin operators that transform in accordance 



with these matrices. The ion spins in one unit cell are 
numbered one to eight as shown in Fig[T] The space in- 
version symmetry, together with the experimental facts 
that Si — S3, S2 = Sa,S§ = Sj, Sq = —Ss and that 
the spin moments are only in a-b plane, suggests that 
the possibil ities of magnetoelectric coupling term can be 




FIG. 1: A sketch of spin structures from a top view along the 
c axis in one unit cell of RMniO^. The red lines reflect the 
frustrated magnetic coupling between two chains. 



narrowed down to 

H em = i\ x P y (-S%(q)SI*(q) + Sf(q)S^(q) ~ cc) 
+ i\ y P y (-S%( q )Sr(q)+SUq)Sr(q)-c.c). (2) 

Converting these spin operators to the two antiferromag- 
netic orders, we can simplify EqJ5]to EqfT] 

Effective magnetic model: In condensed matter 
physics, an effective model at low energy is largely in- 
dependent of microscopic models if they share same es- 
sential physics. Therefore, one can derive the effective 
model on a much simplified lattice structure. In the case 
of RMn 2 Os, the important magnetic physics along the 
a axis are two antiferromagnetic chains with frustrated 
coupling^? ? ]. Experiments have shown a 1/4 commen- 
surate magnetic wavevector along the c axis, which can 
also be viewed as an antiferromagnetic order if the unit 
cell is doubled along the c-axis. Therefore, we can derive 
the effective model from a microscopic model with two 
antiferromagnetic orders defined on the two interpene- 
trating sublattices as illustrated in Fig[21 where </i,(2) are 
the effective antiferromagnetic exchange couplings which 
establish two antiferromagnetic orders, J 3 is the effective 
frustrated coupling between two chains in one unit cell 
along the a axis and J4 is the effective frustrated coupling 
between two chains in two neighbor unit cells along the 
c axis. Using standard field theory [? ? ], we can show 
that the effective field theory described by the two anti- 
ferromagnetic orders n\ and n 2 is given by the following 
Hamiltonian 
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FIG. 2: A sketch of lattice structures of two frustrated coupled 
antiferromagnetic Heisenberg models in two dimension. 



where g{T) = go + g\T is a temperature dependent 
parameter induced by the quantum and thermal fluc- 
tuation, the parameter a reflects that the intra frus- 
trated coupling in one unit cell is larger than the in- 
ter frustrated coupling along the a axis between two 
neighbor unit cells along the c axis, and the param- 
eter D describes a possible magnetic easy axis along 
the a axis. From the microscopic coupling parame- 

2A-, 
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ter, g x = 0.26(Ji + J 2 )Sa~ 2 {^^ )*, g 2 - «.-* (Jl+Ja)s , 

n — 7 C"2 _ T C2 n _ '(Jl+J2)S 2 (J3-J4) „ n j 

pis — Jib , p 2s — J 2 b , a — 1 (j 1+ j 2 ) a ana 

v = {.h+.h)S 2 Adding the lattice dynamics, we 

reach the total effective Hamiltonian as 



2M^ y) 



H„ 



(4) 



where ir y is the conjugate momentum of P y . 

The Hamiltonian in Eqf5] precisely captures the phase 
diagrams of RMn 2 0§[l ? ? ]. To study the mag- 
netic phase diagram, we can integrate out the lattice 
dynamics. After integrating out the lattice dynamics, 
the effective magnetic Hamiltonian is the Hamiltonian in 
Eq 11 01 with a replacement of g by g(T) = go + g\T where 



5o =9o+87, D by D 



D 
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and an additional 



term — ^y(nfn 2 — nfn^) where 7 = g=- and X± = X x ±X y . 
It is clear that the a term favors an ICM phase while 
g(T),j and D favor a CM phase. If D ^ 0, at relatively 
high transition temperature T c ± , the model exhibits a 
first phase transition to a collinear magnetic phase with 
order < nfn 2 >^ and then exhibits a second phase 
transition at T c2 < T c \ with order < rv[n v 2 >^ 0. The 
collinear phase becomes a col-plane phase. The ferroelec- 

tricity is given by < P y >= — — - — 3 — 2— . With 

the proper values of a, at a low temperature Tic, the 
ICM phase can win over the col-plane magnetic phase. 



In the ICM phase, the global average, < P y >— 0. Fig[3] 
sketches the phase diagram. The phase diagram quali- 
tatively matches the current experimental results on the 
phase diagram of RMn 2 0^ ? ? ]• Results from mean 
field or large N limit calculation, using the real experi- 
mental data as input, will be reported elsewhere [? ]. In 
this paper, we present a thorough study of the simplified 
version of the model in one dimension which captures the 
essential low energy physics. 

Coupled Sine- Gordon model: We focus on an one- 
dimensional effective model by ignoring the dynamics in 
z direction in Eqf4] The one dimensional model still cap- 
tures essential physics of the frustration. For the one 
dimensional model, we can go beyond mean field calcula- 
tion and evaluate the dynamics in a controllable pertur- 
bation manner. Using the angle to parameterize the an- 
tiferromagnetic vector in x-y plane as fii = (cos9i, sinOi), 
we obtain the effective Hamiltonian in one dimension as 
a coupled Sine-Gordon model, 



H 
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2^ dx ' 
2 [ dx ■ 



a sin 6L — 

ox 

D cos 9+ cos 0_ 



g cos 



7 cos 



(5) 



where 9± = 9\ ± 9 2 . Minimizing the potential part in 
Eqf5l the CM-ICM phase transition line is estimated to 
be 



p(g(T) +1 /2 + D) = a 2 /2 + 



D 2 



2{c?-pg{T)y 



(6) 



To obtain a more accurate result, we can apply a 
controllable numerical method[? ]. In the CM phase, 
the ground state of the model is always described by 
6* + = 6*_ = or 7T. In the ICM phase, it is clear 
that the spins on the two antiferromagnetic chains are 
both rotating with the same wave vector q through- 
out the system. Therefore, we can take the follow- 
ing ansatz: 6L(ir) = X^o( a ™ cosngx + b n sinnqx) and 
0+(x) — qx + X)^Lo( c « cos nax + d n sinnqx). The coef- 
ficient in this ansatz decreases rapidly as n increases. 
As an example, for a typical set of parameters {p = 
2, a = 0.4, ff = 0.03,7 = 0.01, d = 0.005}, a nu- 
merical solution of the ICM ground state is given by 
9+ = qx + 0.0322 sin(2 9 a;) - 0.000255 sin(4gx) + ... and 
9- = -§ +0.0526 cosmic) +0.0000676 cos(3ga;)-h.., where 
q = 0.199. Therefore, the first two terms in the ansatz 
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FIG. 3: The phase diagram as go vs T. The vertical dashed 
line represents phase transitions in most of KMn^O^, mate- 
rials. The inset (a) is the ferroelectricity as a function of 
temperature. 



can give us a very good approximation of the ground 
state. 

Energy dispersion of magnons in the CM phase: In the 
CM phase, the ferroelectricity, < P y >= — A+ 2^ A ~ ■ The 
dynamics of magnons and the dynamics of electric de- 
gree of freedoms are decoupled. To show this, one can 
expand the free energy at the vicinity of the ground state, 
u\ = 59 '_(_, u 2 = 59-, [3 — 6P y and show that there is no 
second order coupling between Ui and (3 (the lowest order 
coupling between m and (3 is \ + [3u 2 +X-f3uf). Therefore, 
there is no electromagnon excitation. The energy disper- 
sions of the two magnons in the CM phase are given by 



El 



\J p{pk 2 + 7 + .g + D± \/{g-"i) 2 + a 2 k 2 ). (7) 



Electromagnons in the ICM phase: In the ICM phase, 
< P y (x) >= — ^cos(qx) where q is the ICM wavector. 
The magnon and phonon dynamics are coupled due to 
the ICM modulation. Therefore, the low energy excita- 
tions are electromagnons. To calculate the dispersion, we 
expand the Hamiltonian at the vicinity of the ICM phase, 
ui = 59 + , U2 = 59- and (3 = 5P y . The fluctuation up to 
the second order of these dynamic variables is given by 



5H^ = SH { > + 5H\ >. The first term is the free part 
given by 

2 

fc 8=1 

with E 1 (k) = pk 2 +°§±+ A - (ai 4 A K +A - } , E 2 (k) = E 1 (k) + 

A 2 

Ao where Ao = aq — 2g — and the phonon frequency 
— \ffi (for convenience, we have taken the mass M = 1 
in Eq0]) . The second term is the interaction part 

5E[[ 2S> — J dx[aai cos(qx)u\u 2 + D sm(qx)uiu 2 



+~(3{—\-U\sin{qx) + X + u 2 )]. 



(9) 



In EqUJl the terms in the first line couple iti(fc) with 
u 2 (k± q) and vice versa and modify the gap Aq between 



two magnon modes. In general, as shown in Fig[4] the 
coupling creates a distinguished kink in dispersion curve 
around k — q/2 rather than a finite energy jump at 
k = q/2 which is normally expected in the ICM phase 
described in a simple Sine-Gordon model[? ]. The terms 
in the second line describe the coupling between the 
magnon and the phonon. The coupling results in two 
general effects. First, the coupling leads to a change of 
phonon frequency. Under condition Wg >> pEi(0), the 
shift frequency of the phonon is roughly given by 5ui 2 = 

4(^-pB 2 (0)) + S(^-pE 2 (q)) - _ 

us to measure the magnons in optical conductivity. The 
optical conductivity is given by o~(tu) = ujhnlGppiu), 0)], 
where Gpp(ui,k) is the full propagator of the ft. In our 
model, we expect double peaks in the optical conductiv- 
ity at the gap energy of the two magnons in the ICM 
phase. Up to the second order, we have, Im[Gpp(ui, 0)] = 

7r[p§z5(u 2 P E 1 (q)) + p^5{u 2 



Second, the coupling allows 



pE 2 (0)) + (1 - 
In Fig (4) by numeri- 



cally solving the dynamical equations of 5H^ 2 \ we plot 
the result of cr(w) and the dispersions of electromagnons 
with a typical parameter setting {p, a, A+, A- , k, D, g} = 
{2, 0.4, 0.06, 0.06, 2, 0.002, 0.005}. 

Effects of external magnetic field: In our model, since 
n is the staggered moment field, the effect of an external 
magnetic field H is equivalent to creating an easy plane. 
However, since the staggered moment is in the a-b plane, 
the effect of the field is only important when the field is 
along the a or b axis. In the presence of the easy axis 
parameter D along the a axis, the effects of the external 
magnetic fields along the a axis, if a , and b axis, Hi,, have 
exact opposite effects. They simply change D to 



D(H a , H b ) = D- H 2 a /2p + Hl/2p 



(10) 



in Eqd From E x (k), E 2 {k), Eq(9]and EqEl EqllO] leads 
to a few important and immediate predictions. First, the 
effects of H a and Ht do not depend on their directions 
along their own axis. Second, Hb can drive the system 
from the ICM phase to the CM phase while H a can drive 
the system from the CM phase to the ICM phase. This 
result has been observed experimentally in [? ] . The crit- 
ical fields to drive the transition can be estimated for the 
one dimensional system from Eq|SJ Third, in the ICM 
phase, the energy dispersions of the electromagnons as a 
function of H a and iff, can be predicted. From E\(k), 
E 2 (k) and EqEl the energy of the electromagnons is ex- 
pected to increase (decrease) as Hb (H a ) increases. Fi- 
nally, the external magnetic field does not add additional 
peaks, which contradicts the conventional picture of the 
Zeeman energy splitting of magnons. 

In conclusion, we develop an effective model that ex- 
plains the phase diagram and the mechanism of mag- 
netoelectric coupling in multiferroics RMn 2 0^. To our 
knowledge, this is the first theoretical effective model for 
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FIG. 4: The numerical result of the optical con- 
ductivity versus frequency with the parameters set to: 
{p,a, A+, X-,K,D,g} = {2,0.4,0.06,0.06,2,0.002,0.005}. 
The inset is the dispersion of the electromagons. 



these materials. A detailed study of low energy excita- 
tions in one dimension is performed to explain the selec- 
tion rules of electromagnon in optical conductivity mea- 
surements [? ]. Our prediction of the electromagnon dis- 
persion and its dependence on the external magnetic field 
can be tested in future experiments. A quantitative study 
of our model incorporating the experimental data on dif- 
ferent RMn^O^ materials will be reported elsewhere[? ]. 
We expect that the model presented here can be applied 
to other multiferronics materials where ferroelectricity is 
correlated to a collinear or col-plane magnetic phase. 
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